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Abstract 

The 2 — > n scattering with final particles at rest is discussed. The 
comparison with purely soft processes allows to identify symmetries 
responsible for vanishing of certain 2 — * n amplitudes. Some examples 
are given. 
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1 Introduction 



In the present paper, third of the series (cf. [T], 0), we consider 2 — > n 
processes: two particles of opposite momenta scatter to produce n bosons at 
rest. This process has been already considered, on the tree level, using prop- 
agator approach [3] -j- jH], diagrammatics [Sj or Feynman wave function 
method 0, PQ. The last method seems to be the most efficient. 

In Sec. II we present a very simple and straightforward derivation of the 
Brown-Zhai [9 a algorithm. We give a general procedure to construct 2 —>■ n 
amplitudes which, in particular, allows us, for a wide class of theories, to 
compare 2 — > n amplitude with a purely soft one. By using the results 
concerning the relations between symmetries and vanishing of purely soft 
amplitudes jTUj, |2] we are able to identify, in certain cases, the symmetries 
underlying nullification in 2 — > n scattering. 

In Sec. Ill we discuss some examples involving only bosons while Sec. IV 
is devoted to fermion-antifermion scattering. The latter is discussed in more 
detail and the complete expression for the amplitude is given and compared, 
for n = 2, with diagrammatic result. 

Sec. V contains some conclusions. 
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2 Calculating 2 —> n amplitudes 



Let us first remind the general procedure for calculating the matrix ele- 
ments of field operators in the tree approximation j^j, PQ. Our starting point 
is a lagrangian 



1 - 

L = - Y^id^iPQi - m 2 ® 2 ) - V®), (!) 



where {$i}^Li is a multiplet of scalar fields; higher-spin fields can be also 
included. 

Let {fn }™=o be a complete set of normalized positive-energy solutions to 
the Klein-Gordon equation of mass m 8 ; in what follows we shall use momen- 
tum basis, fff\x) = -. 1 , :e~ ipx , p 2 = m 2 . Define the classical free field 

P yj (2tt) 3 2E p 

by 



M*)=YlWnf!PW + fflf»( X )) ( 2 ) 
n 

Consider the set of integral equations for 

d 4 yA Fi .(x-y)-^-^ (3) 

where Ap..(x) = Ap(x; m 2 )5ij is the Feynman propagator. Under some as- 
sumptions eqs. (JBJ) admit unique solution, at least at the formal perturbative 
level. Once eq. (J3J) is solved we obtain $j( x\ clS cl function of fi n and fin\ 

$i(x) = <$> l {x I fi, fi) (4) 

Taking succesive derivatives with respect to fi's and fi's at fi — fi — one 
obtains the matrix elements of the field operator between in- and out 

states in the tree-level approximation PQ; for example, 



d 2 ^{x\ fi.fi) 

dfi n k) dfi§ 



q =< (ml); out | $i(a;) | (nk);in > 



(5) 

tree v 7 



3 



In practice, it is more convenient to work with differential equations in- 
stead of integral ones. Eqs. (12131) imply 



I 13 J) |„-.„= (7) 

The converse is, in general, not true. Indeed, for example /?'s and /?'s enter 
| /3,/3) as arbitrary constants; they can be replaced, without alterning 
© and (J7J), by any coupling constant dependent functions /3(A), /3(A) such 
that /3(0) = /3, /3(0) = (3. The solution to (jUJ), (jZJ) can be made unique 
by adding further constraints following from the tree-graph interpretation. 
Consider $ 4 -theory as an example, 

vm = ^ A 

Topological relations for tree graphs imply that /3, (3 and A enter the matrix 
elements of field operator in the combination \ n /3 k j3 l with k + I = 2n + 1; in 
other words, the contribution of the A "-order is a homogeneous polynomial in 
(3 and f3 of degree 2n+l. This additional condition makes (0), © equivalent 
to ©. 

Let us now apply the above formalism to the problem of 2 — * n scattering. 
Assume we have two initial particles of the same kind, carried, say, by $i, 
with the momenta p and — p\ we want to compute the amplitude for producing 
n particles at rest, the final particles being of different kind than the initial 
ones. 

Differentiating (JHJ), with respect to (3^ we arrive at the following set 
of equations 

2 . d$j(x) d 2 V d$j{x) 



(i) 



= SnfP(x) (9) 



v=o 



Now, 8$J ( ^ generates, in the tree approximation, the matrix elements of z-th 

9 'ftp 

field operator with one " 1" particle carrying the momentum p already present 
in the initial state. The second initial particle is obtained by applying LSZ 
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reduction procedure to the matrix elements of $i(x). Therefore, we have 

only to generate final particles at the threshold and the only parameters we 

—ft) 

need to keep nonvanishing are (3 Q , % = 2, N. This allows us to simplify 
the whole algorithm. First of all note that all matrix elements generated by 
the solutions to © are now x -independent and eqs. (JHJ), (JTj) can be rewritten 

as 

{dl + m*)* i (t | (3,]3) + = (10) 

d$i(t | (3,(3) 

Ht\P,P)W=0= 7=%E ^ (11) 



Denote 



(12) 

again with all (3's and (3 s, except (3 Q , i ^ 1, vanishing: ipi(x \ (3,(3) generate 
matrix elements with one particle, carrying the momentum p, in the initial 
state and arbitrany number of particles at rest in the final one. Therefore, 
due to the translational invariance, 

^\P,P)=Ut\P,P)e 1 ^ (13) 

and eqs. (JHJ), © take form 

(9? + Mt)Ut \0J3) + M .^ j{t) Mt I 13,8) = (14) 

* w '"■*= ^7wm e " Ml ' • (15) 

where Mf = mf + 

Having solved (|14|). (fTB^) one can compute the relevant amplitude from 
LSZ formalism: 



i(27r) 3 5^(p + q) 

y/(2TT) 3 2M 1 



oo 

tMit/o2 , A/ r2 



dte- lMlt (df + Mf)^(t | p, (3) (16) 
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generates all 2 — * n amplitudes. 

Eqs. (f*T3|) (Unj) provide the general solution to our problem. 

Particulary interesting situation emerges if the theory exhibits additional 
symmetry 

$l = -$i, $ j — $ J = $ j, i^l (17) 

Consider first eqs. ©, ©; once we put (3^ = = (3^ both eq. (0) and the 
initial condition (J2J) are invariant under (ftTJ). By uniqueness of the solution 
we find 

I 0,^1^0^=0 ( 18 ) 

This conclusion is obviously equivalent to the statement that only graphs 
with even number of external "l"-lines are allowed. The counterparts of 
fITUJ) . ([TED read now 

{dl + m?)$ 4 (t) + ^ffi = , i = Z, N (19) 

Mt) \v=o= fot imit ,i = 2, ,iV (20) 

Vl (§L) = vm |* 1=0 , (21) 

which is a set of coupled dynamical equations for the reduced system obtained 
by ignoring $i. Let us further consider eqs.fJHJ). We show that 



(i) 



, 1 (22) 



This is, again, a consequence of the fact that only even number of external 
"l"-lines is admitted. The formal proof is based on symmetry (|17j) : from the 
uniqueness of the solution to (jOJ, (|Zj) it follows then that $j(x | /3,/3), i ^ 1, 

are even functions of /jl 1 ^, which implies f!22j) . 

Taking into account the above property we can rewrite (|14jl 



(a t 2 + M 1 2 )^ 1 (t|A/3) + J^ 



-^(t|/3,/3) = (23) 

>t>l=U 



^(t) | y=0 = e-' Mlf (24) 

y /(2n) 3 2M 1 
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We see that in the presence of the symmetry (|17|) the problem gets simplifield: 

first we solve N-l equations for $j(t | /3, f3), i = 2, ,N and then one 

equation flU for ^(t \ f3,]3). 

Let us now make the following important observation. Take the modified 
theory given by the lagrangian 

1 1 - 

L = -{d^id^ - Ml®l) + - Y^iPfi&Qi ~ ™ 2 i®D ~ V(§), (25) 

i=2 

with Mf = m\ and consider the same process 2 — ► n except that now 
all particles, both initial as well as final, are at rest. Repeating the procedure 
outlined above we conclude that the amplitudes for both processes coincide, 
up to irrelevant delta function expressing momentum conservation. This 
result can again be easily understood on the diagrammar level. Symmetry 
(|17|) implies that there is a unique chain of $i-lines connecting both initial 
particles and the three-momentum flows only through lines of this chain. 

The observation made above allows us to explain, at least in some cases, 
the origin of amplitudes nullification in 2 — >• n processes. In fact, for purely 
threshold amplitudes their nullification may result from symmetry properties 
of the reduced hamiltonian system obtained by ignoring space-dependence 
in the original problem jU]. 

3 Some examples 

3.1 Quartic interaction : 0(2)-theory 

Let us start with the model considered in jH]: 

L = \ E( - ™^ 2 ) - ^(*? + $2) (26) 

i=l 

The reduced theory is integrable and separable in elliptic coordinates 

in fact, it is an example of the so-called Gamier system [TI]. In order to 

calculate amplitudes "1" — > "2" at the threshold we look for the solution to 
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the field equations obeying 



a -im\t 

* i(t) i«= imm s zdt) (27) 

a p im 2 t 

i-= Jtm s z ' Jt) (28) 

They read lHj 

$1 = - 7^ 2 2 )(1 - - -^4 + f f ^ i^)" 1 (29) 

V 48ml 48m 2 48m 2 48 2 m 2 m 2 1 27 v 7 

$2 = 22(1 + 7^2^1 2 )(1 - ~ T^~~2 Z 2 + jJr i o ^i^)" 1 ( 30 ) 

v 48mf 17 v 48mf 1 48m 2 , 2 48 2 m 2 m 2 , 1 27 * 7 

where k = 

mi —mi 

In order to calculate 2 — > n amplitude, with all particles at rest, we 
compute 

-a XkzI )(i x ^ e ' imit (K) 

On the other hand consider 2 — > n process with hard initial particles. The 
relevant amplitude can be computed according to the scheme presended in 
[Oj. The relevant wave function obeys eq. (3.3) of [H] with a=2 and is given 
by eq. (3.8) therein 

^(x) = e-^F(l,-l,l-^ ;2/ ) (32) 



A 



48m 2 , 1 - *z? 



(33) 



2 1 48^f^2 



or, explicitly 



A / E p +m 2 \ 2' 
mini \Ep-m2 > 2 



^(a;) = e-*" ( _ 2 7"7 ) (34) 



Eqs. (13 1|) and (|34|) can be now compared. According to the general reasoning 
presented in Sec. II one should identify E p <-> m\ and neglect space-dependent 
factor in (JMJ). Then we see that (}3~T]) and (J53j) indeed coincide (one should 
also take into account the difference in wave function normalization). 

In particular, the nuliffication of 2 — ► n hard-soft amplitudes results from 
Ward identities for purely soft ones [2|. 
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3.2 Henon-Heiles system 

Consider the theory given by the lagrangian 

-§^ 2 " ^ (35) 

It is symetric under $x — > — $ 1; $ 2 — > $ 2 . There exists no stable ground 
state here but on the perturbative level theory is perfectly well defined. We 
shall analyse its threshold behaviour in few steps. 
(A) 2 — > n amplitudes: $i -particles 

We assume that both initial particles are carried by $i field. According 
to the procedure outlined in Sec. II we have to solve first the field equation 
for $ 2 assuming Q 1 = and no space-dependence, 

(d t 2 + m 2 0$ 2 + ^ = (36) 

$ 2 | A=0 = ze im2t = - t e im2t (37) 

The proper solution ( i.e. the one reproducing tree-graph expansion ) is 
obtained by assuming the total energy to vanish; it reads 

ze ini2t 

™ ) = W=&*r m 

This is the generating function for $ 3 -theory. The counterpart of eq. (|2*3*|) 
takes the form 

(D + ml + g$ 2 (t))ip(x) =0 (39) 

Puting 



4j(x) = e-^m, P = < E p ee yjp + m\ (40) 



we arrive at 



{d 2 t - 2z£ p <9 4 + g§2{t))m = (41) 
Ht) \ 9 =o= 1 (42) 
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In order to solve (|41j) we make a change of variable [D] 



x 1 o™^ ° 



which converts (|41|) into 



y(i - v)£ + ((* " — ) - 2y)-f + ^) = o ( i J ) 

dy A m 2 ay A 



The relevant solution reads 



2E 

if>(t) = F(a 1 ,a 2 , 1 p -,y), (45) 

m 2 



where F is the hypergeometric function and 



l± A /l + f* 



«i,2 = ~ 2 (46) 

Let us now apply the LSZ-reduction to produce the second initial particle: 



i 



V(27r) 3 2£, 

i5 {3 \p + q) 
2E p 



d 4 xe- iqx (n + m \)^{x) = 

oo 

J dte- iE *\d 2 t + El){e~ iE ^{t)) (47) 



Now, from (|47)l we see that ip(t) has the following exspansion 

i,(t) = Y,Une mm2t (48) 

n>0 

Taking into account the contribution from n-th term one gets 

= 2m5 {3) (p + g)<5(2£ p - nm 2 ){2E p - nm 2 )a n ; (49) 
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the conribution is nonvanishing only provided a n has a pole at 2E P = nm,2- 
This pole term is easily identified from the solution (J45j) . The result reads 

AmS^(p + q)S(2E p -nm 2 )E p A Tjai + n)T(a 2 + n) 
[ ' ((27r) 3 2m 2 )fn! l 12m| J r( ai )r(a 2 ) 1 ' 

Eq . (|5U|) can be checked by Feynman-graph method. The relevant Feynman 
rules read 



p 2 — vn\ + is 
i 

p 2 — m 2 , + ie 



>9 



-iX 



For n = 1 there is only one contribution coming from the elementary 
vertex 



which contributes an amount IH^E^ m 2 )s < - Hp +g) - m a cr reem ent with (1501) 
for n = 1. 

The case n = 2 is slightly more complicated . There are three graphs 
contributing 



+ + 
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and they sum up to 



(2^J { 9 -6j (51) 

which again agrees with general formula (|50|). 
Let us now note the following: for 

£ = N(N + l) 

A 12 K ' 

one finds a 2 = —N. Therefore, due to T(a 2 + n) /T{a2) = (—N)(l — N)...(n — 
1 — N), all amplitudes A(2 — > n), n > N + 1, vanish. In particular, for 
N — 1 the only nonvanishing amplitude corresponds to n = 1. This latter 
result can be understood by invoking tyhe equivalence with purely threshold 
amplitudes, as discussed in Sec. II. In fact, N = 1 implies, through 
g = ~; in this case the Henon-Heiles system becomes integrable and separable 
in parabolic coordinates [TT]. Let us consider the relevant solution in more 
detail. 

(B) Henon-Heiles model: integrable case 

If g = ^ the reduced dynamics becomes separable in parabolic coordinates 



$2 = -4£ V 



$2 = C + V + « 



^.(Aml-ml) (53) 



The additional integral of motion resulting from the separation of variables 
reads 

F = ^*,(*,*2 - <Mi) + ^(*? + + ^*;*2 + 

(54) 

F generates the following symmetry transformation 

$i -> $i = $i + e\(k$i + ^($i$ 2 - 2$ 2 $i)) (55) 
. £:A 

$ 2 _» $ 2 = $ 2 + 
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We are looking for the solution obeying 



$1 U=o= z l e~ m ^ 



(56) 



v / (2vr) 3 2m i 
$2 U=o= z 2 e m2t 

It corresponds to E = and F = and reads 

$ = 2mi(4m| - m|)(l - apy 

1 A(2m 1 (l-x)(l-i/ 2 )+m 2 (l + x)(l + y 2 )) 1 j 

4/t(m^(l + y 2 ) 2 + 4m 2 y 2 (l-x) 2 ) 

2 (2m 1 (l-x)(l-|/ 2 )+m 2 (l + x)(l + i/ 2 )) 2 1 ' 



x=Zj \( 2m 1 + m 2 ]eim2t 



where 



m 2 \2rrix — m 2 

1/ = \ (60) 

2mi(2mi — m 2 ) 

Let us analyse this solution in some detail. For 1m\ ^ m 2 there are no 
resonances so the relevant threshold amplitudes vanish. On the other hand, 



I / 2m 1 +m 2 \ Z2X im2t , 

a$i _ X ~ {2^T 2 )^ e 2Xz 2 e^ t a 

Ul= °" 1 - SA e ^* " 7^(27^ - m 2 ) + U[Z2) ' ' 

m 2 

which gives nonzero amplitude 2 — > 1 for 2mi = m 2 . This is in agreement 
with general arguments concerning the relation between symmetries and am- 
plitude nullification ^U|, [2|; in fact, for 2mi = m 2 the relevant linear part 
in the transformation formulae ()55|) is absent. 

By replacing mi by E p = \/ p 2 + our expression (jUTjl coincides with 
(|45|) calculated for N — 1. Therefore, the symmetry (J55|) . which is valid for 
arbitrary masses, inplies the vanishing of 2 — > n hard- soft amplitudes. 
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4 Fermions 



It is straightforward to include higher-spin bosons into our scheme. As 
far as fermions are concerned the only but crucial modification is to replace 
the parameters (3, (3 entering (j2J by anticommuting Grassman variables; this 
allows to implement Pauli exlusion principle. 

As an example we shall compute the threshold amplitudes for n-boson 
production by fermion-antifermion pair. The theory we are considering pro- 
vides a toy model for fermionic mass generation is the standard model via 
spontaneous symmetry breaking. It has been considered from the point of 
view of amplitudes nullification in Refs. jl], jS]; diagrammatic approach 
through recurrence relations has been developed in [7j. 

Our model contains massless fermionic field coupled by Yukawa term to 
the Higgs field. The relevant lagrangian reads 

L = iip^d^ + -{d^d**® + m 2 $ 2 ) - - g®^ (62) 
Due to the wrong sign of mass term the vacuum solution is nontrivial 



^ = o, i> = 0, $ = v = V-y- (63) 
Define the physical field p by 

p = $ - v (64) 
In terms of this new field the lagrangian reads 

l = Qiyfdp - M)ij + \{d,pd» P - m y) - ^ P 3 - + 
-gM\ ( 65 ) 

here M = gv } m p are fermionic and bosonic masses, respectively. 

L is invariant under ip — > —ip, ip — > —ip, p — > p; therefore, one can 
apply the strategy of Sec. II. First, we solve field equations for p under the 
condition ip = if) = and 

Q p im p t 

p | A=0 = "7=== = ^ %mpt = z(t) (66) 
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The solution reads 



= TZW < 67) 

1 2d 



or 



<&(*) = „_ J| (68) 

The counterpart of ()23|) takes now the form 

(i-y»d ft -g$(t))il> pa (x) = (69) 
I M 

^)\/=J M =]]j2^E p ^ s >- ipX (70) 
The same reasoning as in bosonic case suggests the following Ansatz for ip ps \ 



r,M = \^Ht))^ ^ (71) 



M 

u{p } s) (72) 



where 



r]J y (2tt)*E p 

is the standard solution to Dirac equation Inserting the Ansatz (|71j) . 
(J72*Jl into eq. (jHEJ) we arrive at 

(£ p + id t - g$(t))F(t)( + F(t) Pk a k V = (73) 
(-£ p - ^ - g$(t))F(t)ri - F(t)p k a k ( = (74) 



or, using (|72|). 



(E p + i$ - g$(t))F(t) — (E — M)F(t) = (75) 
(£„ + zd t + 9$(t))F(t) -{E + M)F(t) = 0; (76) 



the boundary condition (fTUj) is obeyed provided 



| g=0 = 1, F | 9 =o = 1 (77) 

gv — M gv = AI 
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It is not difficult to find the solution to (J75J), ()76)) obeying (|77j). We apply 
(E p + id t + </$(*)) to (J75J) and use (J7BJ); then 

(d 2 - 2iE p d t + g 2 $ 2 (t) + igd t $(t) + M 2 )F = (78) 

Changing the variable, 



2d 



1 _ £W 

1 2« 



(79) 



we arrive at the hipergeometric equation. Its solution reads 



F(t) = F(a, [3, 7 ; y) (80) 

2M , 2M 2E V . , 

a = , /3 = 1 ,7 = 1 " 81 

m p m p vtip 

For F(t) we obtain, respectively, 

F(t)=F(a + l, 7; V)\ (82) 

(JHOI) and JH2J) satisfy (j75 Jl . IfTRjl . 

We are now ready to calculate the exact expression for the tree-graph 

threshold amplitude A(ff — *> np). To this end we reduce the fermionic 
antiparticle by LSZ formula. Then 



- % J d*xv p , s ,(x)(iYd^ - M)ip ps (x) (83) 
becomes the generating functional for all such amplitudes; here 



*M*) = \lj^^V,fy-*' (84) 



In order to calculate (1831) let us write 
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then (J83)) can be rewritten as follows 
-iM 



-iM 
-iM 



d 4 xv{p', s')e- ip x {iYd^ - M)i){x°)e- ipx = 
d 4 xv(p', s')e- i{p+p ' )x {i^p p — M + i 7 °<9oMx c 



E p 



5 (3 \p + p) / dx°e~ 2iE * x v{p', s')(7% — M + i-f°d )<tP(x°) (86) 



Let us note that, due to p ' = —p', 

v(p', s')(7% — M) = 2E p v(p', S ')7° (87) 
and our functional takes the form 

oo 

~ lM 5 {3) (p + tf) [ dx e- 2iEpX °v(p / ,s') 1 (2E p + idoMx ) (88) 



E p 



Consider the last expression. It can be analysed in a similar way as in purely 
bosonic case. In order to isolate the amplitude for the creation of n bosons 
one has to single out the term proportional to Zq . Due to the structure of 
the solution ip ps (x) it can be expanded in Fourier series in e impt . The term 
e inm p t p ro d uces the energy delta function 5{2E P — nm p ); however, because of 
time derivative in (|88J) this delta function is accompanied by 2E p — nm p , so 
the result is zero unless there is an additional pole coming from the gamma 
functions entering hipergeometric series. It is straightforward to isolate this 
term arriving at the following final expression for the amplitude 

A(ff^np)= (89) 
-AmM5{2Ep - nm p )5^(p + f) F(n + 2 ^)T(n - 25£) 



-v(p',s')u(p,s) 



The spin structure of this formula is dictated by invariance properties, angu- 
lar momentum and parity conservation. In order to check the formula (|89j) 
we compare it with Feynman-graph computation. For n = 2 eq. (|89|) takes 
the form 

A(ff - 2p) = 

( M 2 9 M 4 \ 
5{2Ep-2m p )^\p + p)(-^- — )v{p' ) s')u{p ) s) (90) 



Air 2 '' \2m 2 v 2 m p v 2 
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The relevant Feynman rules are 

i( 7 % + M) 
p 2 - M 2 + ie 



p2 _ m 2 _|_ ^ E 



iXv 



Fig.2 

The graphs contributing to A(ff — > 2p) are depicted on Fig. 3 

+ + 
Fig.3 

Again, the amplitude exhibits nullification phenomenon 0], jj], jjj, [5]. 
Due to 

= (1 _ a )(2 - a)...(n - I - a) (91) 

the amplitude vanishes for n > N provided a = — = N is an integer 
(actually, n < N amplitudes also vanish due to energy conservation and the 
properties of spinor wave functions). In particular, for AT = 1 all amplitudes, 
vanish. 
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5 Final remarks 



We presented a simple derivation of the Brown-Zhai Feynman wave func- 
tion algorithm for general theories. It allows to compare our process with 
a purely soft one. The important point is that, for the latter, we can use 
Ward identities following from the integrability of reduced dynamics [2] to 
prove their vanishing. This argument works in purely scalar case (cf. also 
PI) as well as fermionic one ^3]. In this way we identify, in some cases, the 
symmetry underlying nullification of 2 — > n amplitudes. 
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